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We present a calculation of the Ktz transition form factors using the AdS/QCD correspondence. 
We also solidify and extend our ability to calculate quantities in the flavor-broken versions of 
AdS/QCD. The normalization of the form factors is a crucial ingredient for extracting |K, S | from 
data, and the results obtained here agree well with results from chiral perturbation theory and 
lattice gauge theory. The slopes and curvature of the form factors agree well with the data, and 
with what results are available from other methods of calculation. 



I. INTRODUCTION 

In this paper, we consider the extension of the 
anti-de Sitter space/quantum chromodynamics model 
(AdS/QCD) to allow broken flavor symmetry, and ap- 
ply the model to the kaon system and particularly to the 
Ki3 form factors. 

The connection between 5D gravitational theories on 
an anti-de Sitter space and 4D conformal field theories 
began as a correspondence between a type IIB string the- 
ory and an Af = 4 super Yang-Mills theory in the large 
N c limit [EHi. This has inspired an analytic model 
referred to as AdS/QCD connecting 5D theories living 
on an anti-de Sitter space to 4D QCD [J, Q . Interest- 
ing results have been obtained for masses, couplings, and 
electromagnetic and gravitational form factors for vector, 
axial vector and pseudoscalar mesons. For a selection of 
these results, see H 0, 1, S M HO, d EI HI El, M EI 

El, EI EI El, Hi- 

The Ki 3 form factors describe the decays K — > ntv and 
are currently most known for the role they play in high- 
precision extractions of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element V us . The form factors them- 
selves are the K to ir transition matrix elements of the 
strangeness changing vector current. There are two Kt$ 
form factors, bespeaking the fact that the strangeness 
changing current is not conserved, and has a longitudi- 
nal as well as a transverse part. Experiments measure 
the product of the form factors and \V US \. Hence the 
extraction of \y us \ from the AT^-decay measurements de- 
pends on having a reliable calculation of the form factor 
normalization. So far, these calculations have been from 
chiral perturbation th eory l23l [243 . l25l . |2o| or from lattice 
gauge theory [2?], [H, |29| . I30L l3l|. Here we will present 
a first calculation of these form factors from AdS/QCD. 
We will also solidify and extend our ability to calculate 
quantities in the flavor-broken versions of AdS/QCD. 

Our form factor normalizations can be compared to 
those obtained from other calculational methods, our 
slopes can be compared to data as well as to other calcu- 
lations, and since we calculate using an analytic method, 
we can also obtain a curvature that can be compared to 
experimental data. All comparisons of results to other 
methods and to data show good agreement, as will be 
detailed below. 



We work with general mass pseudoscalar mesons. Pre- 
vious results known to us worked in the chiral limit or cal- 
culated some quantities using expansions valid at small 
mass. In particular, while we can neatly derive the Gell- 
Mann-Oakes-Renner (GOR) [13] relation in the chiral 
limit, we do not use it to obtain any of our results and 
can test to see how accurate it is at given quark mass. 

Section HT1 reviews AdS/QCD with notation pertinent 
to several quarks of differing masses. Section Hill gives re- 
sults obtained from two-point functions, including bulk- 
to-boundary propagators, masses, and decay constants, 
focusing on differences from the flavor symmetric case, 
in particular considering scalar mesons along with vec- 
tor, axial vector, and pseudoscalar mesons. Section IIVI 
contains algebraic results for the form factors, with the 
numerical results given in Sec. [V] Closing remarks are 
offered in Sec. IVII 



II. THE ADS/QCD MODEL 

We will use the following metric for the 5 dimensional 
Anti-de Sitter space 

ds 2 = — (ri^dx^ 1 dx v — dz 2 ) , e < z < zq, (1) 



where the metric 
rin V =diag(l,-l,- 



of the 4 dimensional flat space is 
1,-1). The cut-off at z — e (with 
e — > implied) corresponds to UV cut-off in QCD, while 
the hard-wall cut-off at zq corresponds to IR cut-off, 
A-qcd, to simulate confinement. 

According to the AdS/CFT correspondence, for every 
operator in 4 dimensional theory there is a corresponding 
field in the AdSs space. Operators of our interest are cur- 
rent operators J a hji = ip qL 7^ a ip qL , Jj^ = ^q R l^t a ip qR 
and quark bilinear ^q L ^ qR - In the AdSs space, these op- 
erators correspond to gauge fields L®(x, z), RfXx, z), and 
a scalar field X(x, z) respectively. Following jj], we will 
consider a 5D action with SU(3)l<B> SU(3)a symmetry as 
follows 



S = J d 5 x^gTr{\DX\ 2 +3\X\ 2 -^(F { 



F (R) 



The field strength is defined by F i 
i\Lm,Ln\ and analogously for F j 
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and gauge fields interact through the covariant derivative 
DmX = 8mX — ILmX + iXR]\j in such a way that the 
action is chiral invariant. We also use the vector and the 
axial- vector field defined from L = V + A and R = V — A. 

The theory begins as one that has SU(3)l<8> SU(3).r 
symmetry, and one would like to maintain as much of the 
symmetry as possible even when going to massive quarks 
and in particular to flavor non-symmetric quark masses. 
In a chirally symmetric world, the action is invariant as 
the X field transforms via 



X 



(3) 



-> X' = U L XU 
One can expand X as 

X(x, z) = e" a(x ' z)t °X (z)e"" (l ' 2)t ° (4) 

whereupon an axial transformation (which has \j\ = U r) 
induces to leading order a shift in the pion field, 7r' a — 
TT a — 6 a , where 9 a is a parameter in the transformation 
Ul = e~' ie * and is consistent with 7r a being a pseu- 
doscalar field. With flavor symmetry, Xq is a multiple 
of the unit matrix and commuting, so one can easily 
write X = e 2l7r * Xq, as has often been done. How- 
ever, in the flavor non-symmetric case, this would make 
7r a appear to be associated with left-handed transforma- 
tions rather than with axial transformations, and gives 
it unexpected parity properties and mixing with vector 
as well as with axial vector fields. For example, one ob- 
tains a quadratic term in the Lagrangian proportional to 
7 1 MN Ti[X ,d M TT b t b }[X , V N ], which will violate 4D par- 
ity when Xq is not proportional to the unit matrix. (With 
the split exponential, one gets cancellations rather than 
just commutators.) 

Shock and Wu [H| have early on considered three- 
flavor extensions of AdS/QCD, although keeping X = 
e 2l7T * Xq. They did not study the more dynamical quan- 
tities like form factors, but did obtain many good results 
for masses and decay constants. However, as they them- 
selves point out, they did not obtain good results for the 
ground state pseudoscalar kaons with the same parame- 
ters that gave good results for the more excited strange 
mesons. With the exponential split, as we think it should 
be, one obtains good results for pseudoscalar as well as 
strange axial and vect or ( and even scalar) meson states. 

Katz and Schwartz [3J] also considered flavor-broken 
AdS/QCD, although their main focus was on the U(l) 
problem and also did not study form factors. They also 
kept X = e 2 '" 7 "*" ' Xq, but only explicitly studied the part 
of the action that mixes the axial vectors with the pseu- 
doscalars, where problems do not appear. We might re- 
mark already that they used the GOR relation to get 
the strange quark mass. The GOR relation becomes less 
valid as the quark mass increases, and using a different 
method to fix the strange quark mass, we find a larger 
value than they quote. 



Hambyc et al. 



351 ] also studied three-flavor AdS/QCD, 
focusing on quantities that are calculated from four-point 
functions such as the purely hadronic K^i decays or the 



Bk parameter needed to calculate K°-K° mixing. They 
work in a limit where all quarks are massless, and so have 
Xq = X = 0. Hence their subjects and their approxi- 
mations do not greatly overlap with the present paper, 
although we plan to consider quantities obtained from 
four-point functions in future work. 

Turning off all fields except Xq(z) and solving the 
equation of motion, one obtains 



2X 0ij = Vij(z) 



(MijZ + ^Ey; 



(5) 



where ( is a rescaling parameter [36|, l37f discussed below. 
From the AdS/CFT correspondence My can be identified 
as quark mass matrix which responsible for the explicit 
breaking of the chiral symmetry of QCD and Ey as the 
quark condensate (<?i<y) which spontaneously break the 
chiral symmetry of QCD to SU(3)y. Assuming u and d 
symmetry, we have 
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m q 









f °q 










m q 























V o 








(6) 



In general a s =/= a q . However, we will also consider the 
limiting case when a s = a q , as an analytic solution for 
the vector field can be obtained. 

Regarding the quark masses and the quark condensate 
parameter a, we adopt a normalization parameter as ad- 
vocated in [3q . [3?J ; wherein quark masses are multiplied 
by a factor £ = ^/Wcjl-K [38| compared to earlier con- 
ventions and a is divided by the same factor. One can, 
of course, view the quark masses and a in AdS/QCD as 
parameters of this particular model, and many impor- 
tant quantities, including the GOR relation and most of 
the results in this paper, are unchanged by this rescaling. 
However, the rescalcd parameters allow a precise connec- 
tion to the two-point correlation function of the quark 
condensate at small distances, which is known from per- 
turbative QCD, and also leads to a better agreement with 
mass parameters at the hadronic scale and with the quark 
condensate as obtained from methods disconnected from 
AdS/QCD. Quark masses obtained from AdS/QCD had 
generally been strikingly low and the quark condensate 
strikingly high, but following [36|, [37| one can argue that 
the disagreement was a matter of having incommensurate 
definitions. 



III. TWO-POINT FUNCTION 



Up to second order in fields, the action can be written 



S = J d 5 x^J^- z ( d MV^-d N V^ 



3 V M 



2z 



-1 



,_(d M A a N - d N A a M ) 2 + -^-P-(d M 7r a - A a M ) 2 



2z 3 



Ml ( > 



(7) 
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where contraction over t]ml is implicit. The mass com- 
binations come from 



±M v 2 6 ab = -Tr [t\X ] [t b ,X Q ] , 
l M f S ab = Tl{ t- j X }{t b ,X } , 



(8) 



or, 



a 2 



MX 



\ (i>s - v q y 



v s ) 



a = 1,2,3 
a = 4,5,6,7 
a = 8, 

a = 1,2,3 
a = 4,5,6,7 
a = 8, 



(9) 



where 



u g (z) = (m q z + -^-z 3 , 
v s (z) = Qm s z + ^-z 3 . 



For later convenience wc define 



a a (z) 



2 A/fa 2 



2 7l,fa2 



(10) 



(11) 



As shown in Eq. ([5]), the vacuum solution contains both 
explicit and spontaneous symmetry breaking parameters, 
Mij and respectively. The parameters m q and m s in 
the 5D theory are usually considered to be explicit sym- 
metry breaking, and give quark mass terms in the 4D 
theory that are also explicit symmetry breaking. The 
condensate parameters may be considered spontaneous 
symmetry breaking, but in the absence of the quark mass 
parameters (i.e., in the chiral limit where the m q and m s 
go to zero), one expects that the condensate parameters 
are all the same. Hence one may argue that the differ- 
ences in the condensate parameters do arise from explicit 
symmetry breaking. Since the M v functions depend only 
on differences v s — v qi one can say that they would be 
zero if there were only spontaneous symmetry breaking. 
In this limit, the masses of the vector mesons in the same 
octet are degenerate. 

The axial sector of action ([7]) is invariant under gauge 
transformation. 



A. Vector sector 

The vector field satisfies the following equation of mo- 
tion 

r} ML d M f ~ (d L V$ - d N V£)\ + = 0. (13) 



For the transverse part, d ,1 V^ ± (a 



0, one obtains 



d z -d z 

z 



q — a 



(14) 



where q is the Fourier variable conjugate to the 4 dimen- 
sional coordinates, x. 

We shall write the vector field in terms of its boundary 
value at UV multiplying a profile function, or bulk-to- 
boundary propagator, V" x (q, z) = V®±(q)V a (q 2 , z), and 
set V a (q 2 ,e) = 1 (Note that there is no summation over 
the group index a of the profile function). The bound- 
ary value V® (q) acts as the Fourier transform of the 
source of the 4D conserved vector current operator. At 
the IR boundary, we choose Neumann boundary condi- 
tion d z V a (q, zq) — 0. In the = aSij limit, the solution 
can be written in terms of Bessel function 



for q 2 > a a , and 



V V, z) = Qz ( M^I.iQz) + Kl (Qz) 



Io(qz ) 



(15) 



(16) 



for q 2 < a a , where q = \J 1 q 2 — a a and Q = y 'a° — q 2 ■ 
Near the UV boundary, the profile function can be writ- 
ten as 



V(q 2 ,z)^l + ^\n(q 2 z 2 ) + 



(17) 



The longitudinal part of the vector field, = 9 M £ a 



and V z are coupled as follows 

-q 2 d z r(q 2 , z) + a a d z ir a (q 2 , z) = , 



(18) 



1 



Q 



8 Z [ -d z( j ) a ( q ^z)j-^-(r(q 2 ,z)-n a (q 2 ,z))=0, 



(19) 



A a 



A a M — 8m A a , 
7r a - A a . 



(12) 



where we define V" = —d z Tr a and £ a = <fi a — 7r a . The con- 
stancy of a a (z) when Sjj = aS^ simplifies above equa- 
tions into 



Hence, we are free to set A a z = 0. For the vector 
sector, the mass term destroys the gauge freedom for 
a = 4, 5, 6, 7. Hence, we can set V" = only for 
a = 0,1,2,3,8. Wc will show that the non-vanishing 
V z is related to the non-vanishing longitudinal part of 
the vector field. 



o z -d z H 

z z 



ev,z) = o. 



(20) 



This is precisely the equation for the transverse part of 
the vector field. Fixing boundary conditions as 4> a (q, e) = 
and Tt a (q,e) = — 1 on the UV brane and Neumann 
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boundary conditions on the IR brane, one concludes that 
in the limit where a s = a q = a the profile function for 
the longitudinal and the transverse part of the vector field 
are identical, £, a (q 2 , z) = V a (q 2 , z), with a solution given 
by Eq. (fTS"]) . In general, this is not the case, and both 
£ a and V a can be solved numerically. For a = 1,2,3,8, 
longitudinal vector fields are unphysical in the sense that 
they can be eliminated by fixing the gauge, V£ = 0. 

Two-point functions can be calculated from the 
AdS/QCD correspondence by evaluating the action ([7]) 
with the classical solution and taking the functional 
derivative over V® twice. One obtains 



1 

0.5 


-0.5 
-1 



/ ' 
/ / 
/ / 
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N N 0.2 0.4 0.6/ 0,:8 

\ \ ' 

\ \ / i 
\ < i 



i [ e iqx (0 \TJt a {x)Jl b {0)\ 0) = -P^5 ab 9zVa( f ' £) , 

% f e i " 35 (o|rjp(x)j|f(o)|o) = - P ^s ab dzr( f 2,£) , 

(21) 



where Pt^ v 



and Pl v 



q^q v ' jq are 



the transverse and longitudinal projector respectively. 
Comparing this result with the quark bubble diagram 
of QCD, one can fix parameter (75 of the model [i| 



.9 5 = 



12^ 



(22) 



Hadrons correspond to normalizable modes of the 5D 
fields. These modes should vanish sufficiently fast near 
the UV brane such that the action is finite and at IR 
brane satisfy Neumann boundary condition. The eigen- 
value, q 2 = M 2 , is the squared mass of the n-th Kaluza 
Klein mode. We expect vector mesons to be normalizable 
modes of equation (|14p and scalar mesons to be normal- 
izable modes of Eqs. (JTSJ) and p9[) . In the = crSij 
limit, the scalar meson has identical mass with the corre- 
sponding vector meson. However, for a = 1, 2, 3, 8, longi- 
tudinal modes are unphysical. Hence, the lightest scalar 
meson obtained from the longitudinal mode of the vector 
field is a strange meson, fQ. Regarding scalar mesons in 
AdS/QCD, see also Ref. [Hg. 

As a remark, we may note that one could include a 
scalar field explicitly by defining X, similarly to [36|, as 

X = e""*" (X + S) e""*" (23) 

with S = S a t a (ignoring the scalar singlet). In this case, 
after some manipulation, the quadratic terms in the ac- 
tion involving V£, and S b for a, b = 4, 5, 6, 7 become 



5* = 



*■ a 

Mf{z) 
2z 3 

S a S a ' 



-1 
45.5 z 



2 f a68c) ( S b 



2z 5 
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MN 



My 



( ^d N M v 



(24) 



FIG. 1: Plot of ipi (solid red curve), %p2 (dashed blue curve) 
and ip3 (dash-dot green curve), with z in units of No- 



where one can show that the square bracket on the last 
line is zero. Redefining the the vector field, 



(25) 



one obtains a massive vector field and eliminates the 
scalar field, and the action becomes like the one we use 
here. 

When the condensate parameters are all the same, 
wave functions for the vector mesons are given by 



V2zJi(zJM~f - a a ) 



(26) 



with normalization condition, j(dz/z)^ 2 = 1. In par- 
ticular, we obtain an infinite tower of KK rho mesons 
for a = 1, 2, 3, the corresponding tower of K* mesons for 
a = 4, 5, 6, 7, and w° mesons for a = 8. The Neumann 
boundary condition on the IR gives Jq(MPzq) — 0. Iden- 
tifying the lightest mode as the rho meson, we can fix z$ 
parameter of the model. 

Figure Q] shows wave functions of the three lightest 
vector mesons. By our choice of boundary conditions 
hadrons wave functions localized closer to IR brane than 
to UV brane. 

The presence of Jo(qzo) on the denominator in Eq. (|15[) 
indicates the existence of poles for timelike q. More pre- 
cisely, the profile function can be written as a sum over 
mesons poles 



V a (q 2 ,z) 



- ff5 F°V>°(z) 
q 2 - Mf 



(27) 



where F° = \d z ^ n {e) / (g b e)\. From Eq.([2T]). F£ can be 
identified as the decay constant of the n-th KK vector 
meson, 



J a ^(0)\v b (q,X)) = F^(q,X)5 ab 



(28) 
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One can substitute the above expansion of the profile 
function into (|21[) and obtain self energy function as a 
sum over narrow mesons poles. A well known signature 
of large N c QCD, which is intrinsic to the AdS/QCD 
correspondence. 



B. Axial sector 

Many of our derivations for the axial sector resemble 
corresponding derivations for the vector sector. For ex- 
ample, the equation satisfied by the transverse part of 
the axial- vector field is analogous to Eq. (|14p . with a a 
replaced by f3 a {z). 

The profile functions of the longitudinal part of the 
axial-vector field and the w field satisfy the following 
equations 



-q 2 d z (j> a {q 2 , z) + (3 a (z)d z TT a (q 2 ,z) = , 



(29) 



-d z <t> a {q 2 ;Z) 



b a (q 2 ,z)-7T a (q 2 ,z))=0, 
(30) 



where longitudinal part of the axial-vector field denoted 
by j4^|i(ar, z) = d fl <p a (x i z). The boundary conditions 
are <j> a {q 2 ,e) = 0, n a (q 2 ,e) = -1, and d z <f) a (q 2 , z ) = 
d z ir a (q 2 , Zq) = 0. Note that these equations as well as 
boundary conditions are analogous to the longitudinal 
part of the vector field. 

In order to solve the coupled equations, one can com- 
bine Eq. (f29|) and (|3Q[) into a second order differen- 
tial equation, defining y a (q 2 ,z) = d z (j) a (q 2 ,z)/z, to ob- 
tain pjl 



d z 



d z y a (q 2 ,z))+z 



P a (z) 



l)y a (q 2 1 z) = 0. 



(31) 

In this notation, the boundary condition at the IR 
limit is given by y a (q 2 ,zo) = which is nothing but 
d z (j> a {q 2 ,z ) = 0. At the UV boundary, Eq. (gOD and 
boundary conditions of it and <fr, give sd z y a (q 2 7 e)/f3(e) 
1 . Near the UV cut-off the profile function can be written 
as 



y a (q 2 , e) = P a (e) ln(ge) + c 2 {qe) 2 \n{qe) 



(32) 



Notice that although the above solution blows up loga- 
rithmically at UV, the profile functions <fi(q 2 , z) as well as 
7r((? 2 , z) do not, because of a multiplication by z in their 
definition. 

Pseudoscalar hadrons in the axial sector are pions, 
kaons and etas. Explicitly, 7r a field can be written as 
follows 



/ '4 

y/2 



7T t 



V2 



2i_ 

x/6 



7T 

K~ 



K+ \ 
^ / 



(33) 




FIG. 2: Plot of 7i"i (lower solid curve, in red), and (f>1 (up- 
per solid curve, in blue) which correspond to ml = 139.6 
MeV, and 7r5 (dash-dot red curve), and cj>2 (dashed blue curve) 
which correspond to = 1892.3 MeV, for a = 1,2,3, with 
z in units of Zq. The units of <f>n and iv® are zZ 1 . 



The corresponding normalizable modes, denoted by 
y^iz), satisfy Eq. (J2J). Their eigenvalues, q 2 = m® 2 , 
are the squared mass of the corresponding hadrons. As 
in the vector sector, the axial sector allows not a single 
mode but an infinite tower of KK modes. These modes 
satisfy yn( z o) = at the IR boundary, consistent with 
d z ^(z ) = 0,<(*b) = 0, and ed z y a n {e)/ f3 a {e) = at the 
UV brane, consistent with 0°(e) = 7r a (e) = 0. Near the 
UV boundary, the normalizable mode behaves like 



y n {z) — a + a 2 z 2 + 



(34) 



or equivalently <j> n ~ a^z 2 /2 and 7r„ ~ rn^aoz 2 / (2(3 a (e)) . 
The coefficcicnt ao is determined by the orthonormality 
condition 



dz 



P a (z) 



y a n (z)y a m (z) 



(35) 



These normalizable wave functions can be solved numer- 
ically. Using this normalization, the dimension of the 
normalizable modes is different from the dimension of 
the profile function. We use it because it is well behaved 
for the ground state in the chiral limit (despite the 1 /m^ 2 
on the right-hand-side). A plot of ir n and 4> n , for n = 1 
and n = 2, is shown in Fig. [5] 

Let us derive how the profile function, y a {q, z), can be 
written as sum over normalizable modes, y^(z). First, 
we write y a (q,z) = J2 c n{l 2 )yl{ z ) ■ Multiplying the 
left and the right hand side of the equation by z(q 2 — 
m c ^)y m {z) I (3 a (z) , then integrating over z, one obtains 



C(<r)- 



P a {z 



■y a m {z)d z y a {q 2 ,z) 



P a (z) 



y a {q\z)d z y a m {z) 



(36) 



after integration by parts and imposing the equation of 
motion (|3 1 [) . The second term and the upper limit of the 
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first term vanish by the boundary conditions of y a (q,z) 
and Un(z). Hence, ignoring the non-pole terms, the pro- 
file functions can be written as 



y 



°(<zV)=x: 



< 2 y°(e)y°(z) 
q 2 — mf 2 



(37) 



which can be integrated to obtain 



q 2 - ml 2 



n * 
TT a (q 2 z) = ^r- 95 " 1 ™ 2 ' 



q 2 — m^ 2 



(38) 



where /° = -d z ^(e)/(g 5 e). 

Axial current-current correlators are analogous to 
Eq. (HU, 

Vxdb d,A a tf,e) 



i I e i « x (0|TJ^(a;)J^_ L (0)|0> = -Pf5 

J X 



gh 



i / e^(0|TJ^j(x)J^|(0)|0) = -P^S 



g¥ 



(39) 



from which, one can identify as the decay constant, 

(40) 

775.5 MeV gives z Q = (322.5 MeV)- 1 . 



^|| (°) Aq) 



The value Mf 
Parameters m 



in i 



and /f, 



i 

for a 



and <7,j can be determined by fitting 



1,2,3, with the pion's mass and 



pion's decay constant respectively. Given experimental 
data for m w = 139.6 MeV and f v = 92.4 MeV, we 
obtain a q = (V3/(2n)) (328.3 MeV) 3 = (213.7 MeV) 3 
and m q = (2tt/V^) 2.29 MeV = 8.31 MeV. In the 



4,5,6,7, with 



<7 q = <7 S limit, fitting mj, for a 
the kaon's mass, rtiK = 495.7 MeV, gives m s ~ 
(2tt/V3) 51.96 MeV = 188.5 MeV (model A). A global 
fit to fifteen observables (model B) yields zq = 
(328.0 MeV)" 1 , m q = (2tt/V3)2.16 MeV = 7.84 MeV, 
a q = ( v / 3/(2tt)) (312.2 MeV) 3 = (203.2 MeV) 3 , 
m s = (2tt/x/3) 56.81 MeV = 206.1 MeV, and a s = 
(\/3/(27r)) (322.8 MeV) 3 = (210.1 MeV) 3 . The fifteen 
obervables include eleven observable in Table [J and the 
additional four observables are /+(0), A^_, A" and Ao, 
from the kaon to pion transition form factor discussed in 
SccJV] The quark masses given here include the normal- 
ization parameter suggested in [H, [13] ■ The AdS/QCD 
quark masses are renormalization scale independent. In 
QCD the quark masses do evolve with renormalization 
scale. We should compare our masses with experimental 
QCD values of the quark masses at a low renormaliza- 
tion scale, say 1 GeV or perhaps a bit below. The quark 
masses quoted by the particle data group [3!| evolved to 
1 GeV using their prescriptions are in the range 3.4-7 
MeV for m q and 95-175 MeV for m s . Predictions of the 
model for masses and decay constants using terms up 



TABLE I: Masses and decay constants. Model A is a four 
parameter fit to four observables as indicated in the Table, 
and maintains a s = cr q . Model B is a five parameter fit to 15 
observables (11 from this Table and 4 from the kaon to pion 
transition form factors discussed in the next section) with 
a s 7^ o q . The values of the parameters are given in the text. 



( 1 ncririfia rilr* 
VJUbcl VdUlc 


Mnrlpl A 

(rr — rr \ 
{(Js — CTq) 


Mod pi R 

(rr -4- rr \ 

(a s 7^ O q ) 


A /To d en fori 




^ivie v ) 


^ivie v ) 


y ivie v ) 




(fit) 


134.3 


139.6 


U 


(at) 


86.6 


92.4 


m K 


(at) 


513.8 


495.7 


Ik 


104 


101 


113 ± 1.4 


m K * 


7Q1 


(SQ7 




fn 


28. 


36 




m p 


(fit) 


788.8 


775.5 


F l/2 


329 


335 


345 ±8 


m K - 


791 


821 


893.8 


F l/2 


329 


337 




m ai 


1366 


1267 


1230 ± 40 


pi/2 


489 


453 


433 ± 13 


m Kl 


1458 


1402 


1272 ± 7 


p l/2 


511 


488 





to second order expansion in the fields of the 5D action 
are summarized in Table HI These may be compared to 
results in [H HI S3] . 



C. Massless pion limit 

The AdS/QCD model has consequences of chiral sym- 
metries, such as the Gell-Mann-Oakes-Renner relation 
(GOR), as shown in [|. Here, wc will present a slightly 
different derivation, starting from the normalization con- 
dition, Eq.([35]). As noted in [|, the weight function 
z/(3(z) has significant support only for z close to z c = 
■\/m q /3a, hence, the normalizable wave function y n can 
be evaluated at z ~ e and moved outside the integral. 
Away from z c , the weight function decays very fast, hence 
the upper limit integral can be replaced by infinity. Not- 
ing that from just after Eq. (|3"8")) . y„(e) = — ffs/jr, wc 
obtain 



glffmV 



dz 



P a {z) 



1. 



(41) 



2m q a q 



The GOR relations immediately follow, f^rn^ 
for a = 1,2,3 and f^m 2 ^ = i m q + ™>s)(&q + c s )/2 for 
a = 4,5,6,7. However, for the kaon case, our results 
deviate by over 30 percent from the GOR relation. 

As approaches zero, the GOR relation becomes 
exact. Fixing /„. to experimental data, parameter a q 
approaches (331.6 MeV) 3 . In this limit, the iri(z) nor- 
malizable wave function becomes constant, tti(z) — 
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— 1/(35/77), throughout the region of interest with a step 
function-like jump near the UV boundary. The wave 
function of the lightest mode can be solved analytically 
in terms of modifed Bessel function. Defining 77 = g^a /3, 
one obtains 



yi = Nz 2 

where, 



7 - M + T7(^iMJ' (42) 



N 2 



^r(2/3)r(i/3) T 



h (v4) 



(v4) 



(43) 



Evaluating y\ at UV boundary, we obtain an equation 
relating and a in the chiral limit which is in perfect 
agreement with previous result (T3 |. One should notice 
that n > 1 KK pions still present in the chiral limit. 



IV. THREE POINT FUNCTIONS AND FORM 
FACTORS 



Matrix elements of the current are obtained from the 
three-point functions using [l2l. [l9|. 

(o)K(fc)) = 



lim 



p — >m ri 

fc 2 — >m? 2 



g^fcg (p 2 -m t ^)(fc 2 -m^) 
p 2 /c 2 



x / d 4 xd 4 we tpx - lkw (0\TJ2i\(x)jf (0)jft(u;)|0) 



(48) 



from which we obtain the form factor /+. 

A similar expression for the longitudinal part of the 
current J/j* (0) using V^(j/), allowing us to obtain the 
scalar form factor /q. 

The relevant part of the action receives contributions 
both from the gauge terms and the chiral terms. Keeping 
only terms that have one vector field and two pion fields 
(cither 4> a (x, z) or ir a (x, z)) one obtains 



The transition form factors for decay is defined 
from [ilEU 

(n-(p)\J^=V\K°(k)) 

= f+(q 2 )(k + P )» + f-(i 2 )(k-p)» , (44) 

with q = k — p. By isospin, they could equally well be 
defined using the ir° transition. Only the vector 

part of the current contributes. Further let 



Mq 2 ) = U(q 2 ) + 



q 



-f-(q 2 ), 



(45) 



L K 



so that /+ and fo come from the transverse and longi- 
tudinal parts, respectively, of jj} AS ^~ 1 \ Unless f-(q 2 ) 
diverges as q 2 — > 0, one has /+(0) = /o(0). One may also 
write 



rf|AS|=l) = Jf + Uf 



(46) 



to show the SU(3) flavor indices. 

In AdS/CFT or AdS/QCD, the three-point functions 
involving three currents are obtained by functionally dif- 
ferentiating the 5D action with respect to their sources, 
which are taken to be boundary values of the 5D fields 
that have the correct quantum numbers 0, d, E3, Eil ■ To 
wit, 

(i/i 3 ) 5S(Virir) 



{o\Tj^( X )jf(y)j^(w)\o) = u ; v ; 

11 x A " 5A^l(x) 5V^(y) 5Afa(w) 

(47) 

where S(Vtttt) is the relevant part of the 5D action eval- 
uated using classical fields that solve the equations of 
motion with z = boundary values AjJ^(x) or Uf^(j/). 



S{Vttit) 



d 5 3 



f" 



be 



x {d»4> a V* v d v (t> c + 2d z d„(f> a V z b d u (f> c ) 

g abc _ 0^o) yb^c ^ 

1 



h abc -d" (7r a 7r c ) - d»<l> a n c V, 



5 ahc <9 z 7r a U z V - h abc -d z (7r a 7r c ) V z b 



The f abc terms come from the gauge part of the orig- 
inal action, and the other terms come from the chiral 
part. We have defined 

g abc = -2iTr{t a ,X }[t b ,{t c ,X }] 

h abc = _ 2l Tr Xq ] {t a^ {tC> X()}} (5Q) 

If none of a, b, or c is equal to "8" , these become 



g abc = jabc^ 



h abc = jabc^ _ Va)Vc 



(51) 



where for Xq = ^cq + cgi 8 



V a = CO + C 8 d a 



a = 1,2,3 



2 (v q + v s ) , a = 4, 5,6, 7 . 



The derivatives indicated in Eq. (|4T|) are facilitated 
by going to Fourier transform space and using the rela- 
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tions hj, \2£ 



<j> a {p,z) 



V z b (q,z) 



n a (j>,z)=7r a (p 2 ,z)^A^ a {p) 



V^(q,z)=V<>(<f,z)V^(q), 

V [ l(q,z) = (4> b (q 2 ,z)-K»(q 2 ,z)) V^(q), 



~d z n b (q 2 ,z) l -^V^{q) 



— > —i (relevant momentum)' 



(53) 



where <\> a and 7r a are now ground state normalizablc 
eigenmodes, with subscript "1" tacit. The superscripts 
on cf> a , 7r a , and V b are the flavor indices for quantities 
with pion or kaon quantum numbers. We are working in 
the isospin conserving limit, so that the 4> a are the same 
for a = 1,2,3 and again the same for the set a = 4, 5, 6, 7, 
and similarly for ir a and V 6 . 



As a check, in the equal mass limit, v s 



the 



transverse form factor should by SU(3) symmetry be 
the same as the electromagnetic form factor. One obtains 
in this limit 



With experience, one can use the above translation dic- 
tionary to obtain form factor results quite quickly. Inci- 
dentally, in the limit of having the same quark conden- 
sate parameter a for all flavors of quarks, one can show 
that the bulk-to-boundary propagator V b (q 2 ,z) for the 
transverse case is identical to (j> b (q 2 , z) — TT b (q 2 , z). 

For the transverse form factor, the V z terms in the 
action, Eq. |49|) . do not contribute. One obtains, 



U(q 2 ) 



dzV 4 (q 2 ,z)\-d z <f> 1 (z)d z ^(z) 



(54) 



0§ 



U(q 2 ) 



dzV{q 2 ,z) 



(55) 



which indeed is the same as F w (g 2 ) as found in Eq. (3.5) 
in [l6[ or to Eq. (38) in [l4j|, allowing for the fact that 
those authors wrote the results using the profile func- 
tions and the massless pion limit, whereas we used the 
normalizable eigensolutions [l9| and nonzero mass. 

The longitudinal form factor is 



f (q 2 ) = jT dz((4>\q 2 , z) - z)) | \d z <fd z <f + g5%( y^V - tt W - J) + |W 



9h 2 



8z 3 (m 2 K - ml) 

<9 2 7t 4 (q 2 ,z) fm'l+ml-ci 2 
(m 2 K - m 2 ,) ' 



(V S - V q )(3v q + Vs)^ 1 - TT 1 )^ 7 - 7T 7 ) - AVqVsCj) 1 ^ 7 - 7T 7 ) + (v q + V s )(3v q - V,)^ 1 - 7T 



(d^ 1 <j> 7 - 1 a z 7 ) 



ml d z a 4 (z) x 7 



8z 3 



2z 3 



2z a 4 (z) 




(56) 



The identity /o(0) = /+(0) is apparent after noting that 
d z Tt b (0, z) = and considering the q 2 = normalizations 
of the profile functions. 



V. NUMERICAL RESULTS FOR K l3 FORM 
FACTORS 

We obtain numerical solutions for the bulk-to- 
boundary propagators and the normalizable cigenfunc- 
tions in the massive quark case using Mathcmatica or 
Maple, and then use the numerical solutions to obtain 



f+( q 2 ) and fo(q 2 )- We present a plot of the results in 

t. [3] Of particular interest for comparison to data (3|| 
and to chiral perturbation theory [H, [H, [H, or 
lattice gauge theory [2?], [H, [29|, [3(| HH are the values of 
/+(0), the slopes of /+ and /o, and the curvature of /+. 

Our results from models A (er, = a q ) and B (cr s inde- 
pendently fit), as well as the results from lattice gauge 
theory, chiral perturbation theory, and experiment are 
listed in Table HU 

Experiments measure /+(0) times the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element V us . If the 
CKM matrix element is gotten from elsewhere, for exam- 
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TABLE II: Results from our models, compared to lattice gauge theory, chiral perturbation theory, and experimental data. 



Observable 


Model A 


Model B 


Lattice 


X PT 


Data [411 


/+(0) 


0.965 


0.936 


0.968(11) [28] 
0.9742(41) [30] 
0.9560(84) [31] 


0.961(8) [23] 
0.978(10) [24] 
0.984(12) [25] 
0.974(11) [26] 




K 


0.0249 


0.0227 


0.0237(23)(21) [31] 




0.0249(11) 


K 


0.0021 


0.0016 






0.0016(5) 


Ac 


0.0123 


0.0140 


0.0128(22)(45) [31] 




0.0134(12) 




0.95 L ■ ■ ■ ■ ■ ■ — 1 

0.02 0.04 0.06 0.08 0.1 0.12 

q 2 (GeV 2 ) 

FIG. 3: The Kiz form factors /+ (solid red line) and fo 
(dashed blue line) plotted vs. q 2 over the physical range per- 
tinent to K — * iiev. The plot is based on the "Model A" 
parameters, where a s = a q . 

pie from the unitarity relation, then all the above values 
of /_)_ (0) are in agreement with experimental data. More 
usually, the calculations are taken to be accurate within 
the stated limits, and are used to extract the most precise 
available values of \V US \ from the data. 

Experiments also measure the slope and curvature of 
the K ft form factors. For /+, both the slope and curva- 
ture can be fit, and are parameterized as [39j 

Ml 2 ) = /+(0) (l + X' + (q 2 /ml + ) + \\>{ ( 9 2 /< + ) 2 ) , 

(57) 

while for /o(<? 2 ) there is a linear fit 

/o(g 2 Wo(0)(l + A (g7"4 + )) . (58) 

Values for the parameters are given in the Table. For the 
experimental data in Table HI1 we took the numbers from 
the FlaviaNet Working Group on Kaon Decays [4l| . 

Additionally, [Htj quotes a result /_(0) = -0.113(12). 
The intercept /_ (0) can be related to the slope parame- 



ters, 

/_(0) = m ^ m ' / + (0) (A -A' + ), (59) 

which leads to /_(0) = -0.141 for model A and /_(0) 

(j.l 10 for model B obtained here or /_(0) = -0.129(18) 
using the FlaviaNet fits to experimental data for Aq and 



VI. CONCLUSION 

We have extended the AdS/QCD model of Ref. @, Q 
to SU(3)l x SU(3)r model with a broken flavor symme- 
try. In order to introduce quarks with differing masses, 
we write the X field with the exponentials of the pseu- 
doscalar field split symmetrically about the classical ex- 
pectation value X . We find that neither the longitudi- 
nal part of the vector field nor V z can be gauged away 
for a = 4, 5, 6, 7. If instead of using the symmetric form, 
one expands X = exp(i27r a f a )Ao, as can certainly be 
done when Xq is a multiple of the identity, the longitu- 
dinal part of the vector field and V z will in general mix 
with the 7T field, hence with the longitudinal part of the 
axial- vector field, and give a parity violating term in the 
Lagrangian. 

We have done both a four parameter and a five param- 
eter fit. The four parameter fit constrains the condensate 
parameter to be flavor symmetric, a s = a q , and the other 
four parameters are fit to the rho meson's mass, pion's 
mass, pion's decay constant and kaon's mass. Predictions 
of the model for the non-dynamical properties of mesons 
such as masses and decay constants are within 20% of 
the experimental data. There is an infinite KK tower of 
pions present, just as there is for vector and axial vector 
mesons. For three-point functions, we have calculated 
kaon-to-pion transition form factors, /+ and fo and ob- 
tain excellent agreement with experiment for the slope 
as well as for the curvature. We further find that the 
intercept, /+(0), agrees very well with lattice and chiral 
perturbation theory calculations. 

The five parameter fit allows a s to vary from a q , and we 
performed a global fit to fifteen observables, including the 
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intercept, slope, and curvature of the K to n transition 
form factors. The results were again good, comparable 
to though somewhat improved as expected compared to 
the four parameter fit. The best value of the strange 
condensate parameter was close to the value in the non- 
strange sector. 

We could perhaps add that we found the intercept 
/+(0) to be somewhat touchy. A drift of either m s or 
<j s away from the best values could lead to a significant 
decrease in its value. On the technical side, in the chiral 
limit, /+(0) becomes normalized to unity. This is because 
the profile functions V b (q 2 , z) in the chiral limit are unity 
at q 2 = for all z and all b, so that /+(0) becomes just 
a wave function normalization integral. However, when 
differing quark masses and differing condensate param- 
eters are used, the profile function at q 2 = is unity 
only for z — and can drift quite far from unity as z 
approaches the IR cutoff, particularly if a s gets far from 
<J q . 

A lack within the AdS/QCD framework is the absence 
of error estimates. On the other hand, in the present 
context, it should be remembered that the form factor 
at q 2 = is fixed by a normalization requirement in the 



equal quark mass limit, so what is really being calculated 
for /+.o(0) is the difference away from unity. Here a 10- 
20% error suffices match the error estimates quoted from 
other methods. 

Possible extensions of present work include calcula- 
tions of four-point functions to obtain the Bk param- 
eter and K-k2 decay amplitudes [HJ, to consider isospin 
breaking in the context of the present model, and to use 
AdS / QCD as a method to study how differing quark mass 
and hence differing pion mass affects the calculated re- 
sults, and compare the trends that are found to lattice 
gauge theory. We hope to return to these topics in future 
work. 
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